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Abstract 

We have calculated the band-/ hybridizations for Ce^^Lai-^^Ms compounds 
(x = 1 and X — > 0; M=Pb, In, Sn, Pd) within the local density approximation 
and fed this into a non-crossing approximation for the Anderson impurity 
model applied to both dilute and concentrated limits. Our calculations pro- 
duce crystalline electric field splittings and Kondo temperatures with trends 
in good agreement with experiment and demonstrate the need for detailed 
electronic structure information on hybridization to describe the diverse be- 
haviors of these Ce compounds. 
75.30.Mb, 71.27.-ha, 75.10.Dg 



1 



Typeset using REVT^ 



A pressing issue in the understanding of strongly interacting electronic materials is how 
to produce realistic theoretical descriptions which encompass both crystalline environment 
and symmetry effects, well treated by ab initio electronic structure theory, with dynami- 
cal effects best treated within many body formalisms. A case in point is heavy fermion 
materials with strongly interacting /-electron states that give rise to huge electronic mass 
enhancements. Some understanding of these systems has been reached in Anderson model 
approaches which assume a nearly atomic limit picture for /-states that hybridize with 
extended conduction states through matrix elements determined from electronic structure 
(local density approximation or LDA) calculations In particular, Gunnarsson and 

Schohammer [|l|, have calculated high energy spectra for a number of cerium based metals 
with a T = variational method, which however ignored crystal field effects. The LDA 
has been used to estimate hybridization induced crystalline electric field (CEF) splittings 
in cerium systems without including the strong correlation effects which give rise to the 
Kondo effect screening of the /-electron magnetic moments by itinerant electrons. Second 
order perturbation theory in the direct Coulomb interaction strength U has been employed, 
and this gives good estimates for electron mass enhancements while only partially capturing 
Kondo effect physics 

In this work, we present first results for a method which combines a nonperturbative, 
finite temperature diagrammatic approach for the Anderson model (the Non-Crossing Ap- 
proximation) with input parameters from the LDA. The NCA can properly generate hy- 
bridization induced CEF splittings while giving an excellent description of the Kondo effect. 
We report calculations of the CEF splittings and Kondo scales in CeM^ (M=Pb, In, 
Sn, Pd) compounds in which experimental Tk values vary with M by nearly three orders 
of magnitude. We have computed energy dependent hybridization matrix elements between 
Ce /-states and other conduction states within the LDA in two limits: 

1) For the dilute alloy system Ce^Lai-xMs with x 0; 

2) For concentrated CeMs compounds. 

We explain and compare the methods and results of those approaches. We correctly find a 
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stable doublet CEF ground state with small Tk values for CePbs and Celna and large 
Tk values with negligible CEF effects for CeSua and CePds. 

In our work, the CEF splittings are induced by band-/ hybridization, which is anticipated 
to be the dominant contribution This splitting arises as follows: a CEF state in Ce- 

configuration is shifted downward by level repulsion through virtual P ^ 

and ^ P ^ P charge fluctuations. In the presence of crystalline anisotropy, different 
irreducible representations (irreps) of the point group in the manifold receive different 
shifts. Using this idea in second order perturbation theory. Wills and Cooper estimated 
the contribution of the hybridization-induced CEF splittings on top of extrapolated point 
charge contributions [0. Although the electrostatic potential from the cubic environment 
can induce the CEF {i.e. , in the point charge model 0), it is difficult to produce a good 
estimate of this contribution in a metal due to conduction electron charge screening [^] and 
metallic covalency. We shall neglect point charge contributions in this letter. 

We describe the Ce^jLai-j^Ms systems in terms of effective impurity Anderson models 
T^,[TB| in the dilute (x — 0) and concentrated (x = 1) limits at a site of cubic symmetry 



relevant to the CusAu structure. In this Letter, we ignore intersite interaction effects such 
as the anti-ferromagnetism found in Celns ||TB[ and CePbs |T^. The impurity Anderson 



Hamiltonian of interest reads 



k(T V s l^cr m 

+ J2^fmfLfm + U ^ UfmUf^., (l) 
™ m<m' 

where m is the label of cubic irrep states, e.g. m = \ J r7+), Vko-m the hopping matrix 
element between conduction electron (ck<j) and /-orbital (/m), Ns the number of sites, e/m 
the /-level energy measured from the Fermi level, U the on-site Coulomb repulsion between 
two /-electrons. 

We have used efm = —2.0 eV for Hund's ground multiplet (J=5/2) in all the calculations, 
consistent with experimental 0] and theoretical [g,^ values. The spin-orbit (SO) splitting 
Aso was read off from the separation between J = 5/2 and J = 7/2 peaks in the Ce 



4/ projected densities of states (DOS). We find Aso = 0.35 eV for all M, in agreement 
with atomic values. We set the onsite Coulomb repulsion f/ ^ oo in our many body 
approximation, though we partially correct for this as we shall describe below. In the cubic 
point group symmetry of the CeMs compounds, the ^ = | multiplet decomposes into a 
Tj magnetic doublet (||;r7)) and Fg quartet (||;r8)), split by an energy A78. The J = ^ 
multiplet splits into Fg and Tj doublets, and a Fg quartet. Experimentally, the F7 doublet 
lies lowest for M=Pb, In (Ayg > 0), while no CEF splitting is resolved for M=Sn, Pd. 

The hybridization matrix elements are calculated from the LDA using the Linearized 
Muffin-Tin orbital (LMTO) method in the Atomic Sphere Approximation (ASA) including 
the so called combined correction term |]10[. We assumed the same Wigner-Seitz radii for Ce 



and M (=Pb, In, Sn) and used experimental lattice constants. For CePds, the Wigner-Seitz 



radius of Ce was set to be 10 % larger than that of Pd |jTT[. We used 165 k-points in the 



irreducible Brillouin zone for the self-consistent solution with DOS integrations carried out 



using the tetrahedron method |jl3[ . We set the orbital basis as s, p, d, f for Ce and s, p, d for 
M-ligands. 

In the concentrated limit (x — > 1, i.e., CeMs), we define the hybridization F™^, (e), in 
terms of an effective impurity Anderson model, with the hybridization derived from the 
overlap between /-orbital at the origin and the rest of the lattice. (See FIG. 1 (a) and FIG. 
2 (a)) More specifically 

TZ^,ie) = -Im Y: VnmV^,^,G'iR, R', e + zr/), (2) 

R,R' 

where V-Rm is the hopping matrix element between m-th /-orbital and ligand orbital at R 
and G'(Il, R', e + irj) the Greens' function of ligand electron created at R and recovered at 
R', with the central /-orbital excluded, as shown in FIG. 1 (a). Now the array of ligand- 
and origin excluded /-orbitals (with origin excluded) serves as an effective static medium 
coupled to the /-orbital at the origin. Our method follows Gunnarsson et. a/.'s suggestion 
to interpret the /-projected DOS as the spectral function of an effective non-interacting 
resonant level model. This corresponds to the first iteration of a "dynamical mean field 



theory" or "local approximation" to the interacting problem ||T5[, which becomes exact in 
infinite spatial dimensions. We obtain the hybridization through Hilbert transformation as 
follows: 
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where G^^f (e ± irj) is the 14x14 matrix Green's function, Pmm^(e) the /-PDOS derived 
from the band calculation with the LDA. We display the calculated x ^ 1 limit model 
hybridization in FIG. 2 (a). 

In the impurity limit (x — > 0), we calculate the hybridization between / and hare ligands, 
that is, 

rZ'(e) = vr E VC„jr^.m'5{t - ek.). (4) 

While, in the concentrated limit, the /-electron hops into already-formed bonding states 
of /(origin excluded)-ligand, this impurity limit hybridization accounts for the overlap of / 
and pre-bonding ligand states. (See FIG. 1 (b) and FIG. 2 (b)) Although this hybridiza- 
tion is calculated for a lattice Anderson Hamiltonian, we can still apply it to impurity 
limit of Ce^Lai-xMs since the Ce-Ce hopping is negligibly small [Q. The procedure for 
computing l^ko-m is to set up a resonant level lattice model Hamiltonian matrix from the 
highly non-orthogonal basis of linear muffin-tin orbitals. We reconstruct the Hamiltonian 



from eigenvectors and eigenvalues of the LMTO equation in the one-center expression ||T0 
Hence, in the one-center basis, the angular momentum states become trivially orthogonal. 
After diagonalizing the ligand sector of the matrix and applying a suitable unitary transfor- 
mation to the new ligand basis, we can directly read off hybridization matrix elements, Vko-m- 
To compute F™^, (e), the Brillouin Zone sum of Eq. (^) is calculated with the tetrahedra 
method [0. In this limit, we used lattice constants of LaMs compounds and readjusted 
the Fermi energy such that the ligand bands are filled up with A^totai ~ 1 electrons (which 
accounts for the lattice- wide removal of the single Ce 4/^ electron). 



As seen in FIG. 2, the x — and x — ^ 1 limit hybridizations are almost identical in the 
high energy region (|e — Ef\ > 0.5 eV). In the low energy region (inset), the peak at —0.1 
eV in the impurity limit ( FIG. 2 (a)) is pushed down to —0.2 eV in the effective medium 
hybridization due to the bonding of Ce-f and Ph-p orbitals. Since the Kondo temperatures 
(Tft-) depend upon the hybridization weight below the Fermi energy (Ep), this bonding 
effect can lead to a completely different scale of Tk as x changes. For CePba [Jl^ and CeSna 



Tsj experimental Tk values are constant with x, and our calculations (Table 1) show this 



within reasonable accuracy, given the exponential sensitivity of Tk to model parameters [|19 
Another difference between the x — > and x — > 1 limits is the extra structure appearing 
in the x ^ 1 calculation above Ep (dashed line of FIG. 1(a)) which is due to flat /-bands. 
Although this feature above Ep is qualitatively different from the x ^ limit, it contributes 
little to A78. 

We solved the U ^ 00 Anderson model by using the well known Non-Crossing Approx- 
imation (NCA) which gives a good quantitative description of Ce compounds except for 
T <Tp <^ Tk (where Tp is a "pathology scale" signalling breakdown of the approximation) 
0. To the flrst order expansion in 1/Ng, with A^^^ the ground multiplet degeneracy, the 
spectral functions of the /° and states are solved for from coupled self consistent non- 
linear integral equations for the p , self energies. To partially correct for our U ^ 00 
approximation, we have estimated the contribution to Ayg arising from virtual occupancy 
between ||; and ||; Fg) by employing second order perturbation theory with U = 6 eV. 
We then added the resulting shifts to the bare /-level positions in the NCA. Ayg values were 
inferred as the peak position separations between ||;r7) and ||;r8) /^ spectral functions, 
which include contributions to all orders in V"^, as per Levy and Zhang and in contrast to 
Wills and Cooper [Q. The Kondo temperatures Tk were interpreted as the low temperature 
splitting between the spectral peak and the peak of the lowest /^ CEF state spectrum. 
Tk roughly depends upon the hybridization strengths as 



JVe: 



where -De// is the effective band width, Tg hybridization strength at Fermi energy for ground 
multiplet, and Ng^N^x) the degeneracy of ground (excited) multiplets. 

The CeMs compounds (M=Pb, In, Sn) have the doublet (||;r7)) states as the lowest 
lying multiplet with Ayg values in good agreement with the experiments [p!6|-p!8|, as listed 
in Table 1. For the heavy fermion systems Celns and CePbs, A78 > (Fy is stable) and 
Ayg ^ Tk, both in agreement with experiment. The CEF splitting comes from the larger 
Fy-hybridization in the energy region from 1.0 to 3.0 eV, which results in the doublet (F7) 
ground state for M=In,Pb. Naively, since Ayg ~ V'^, one would expect the larger crystal 
field splitting for Celns to correspond to a much larger Tk value than for CePbs which 
is not seen experimentally. This common reasoning assumes, however, energy independent 
hybridization, clearly not the case here as shown in FIG. 2. In detail, the hybridization 
strength averaged over energy of Celn^ exceeds that of CePbs which dictates Aygl/^ > 
AtsIp^, while the smaller Celna hybridization near Ep yields a smaller Tk value. Apparently, 
detailed hybridization calculations are critical for quantititative understanding of real heavy 
fermion materials. 

The intermediate valence materials, CeSns and CePds, have A78 <C Tk due to their large 
hybridization, as shown in the table. This agrees with experiment which fails to resolve CEF 
peaks. For CePda, the LDA yields an anomalously huge hybridization (up to 1.5 eV) between 
Pd-d and Ce-f orbitals below the Fermi energy, and as a result, the estimated Tk was nearly 



an order of magnitude larger than the experimental value ||21 |. 

The magnetic susceptibility x{T) provides a measure of the degree of screening of the 
local moments by conduction electrons, for which we show our in FIG. 3(a). Celus and 
CePbs showed CEF effects at about 220 and 50 K, respectively, without sizable moment 
screening until the lowest accessible temperature. Starting with well-localized moment at 
high temperatures, x{T) deviated from Curie- Weiss behavior; x{T) of Celus crossed over 
from ||;r7 + Fg) to ||;F7) magnetic moment regime regaining inverse-T behavior at tem- 
peratures (10~100 K) well below the CEF splitting. For CePbs and Celus, the Kondo 
resonance peak begins to emerge (not shown here) near few K where the numerics of this 



calculation becomes unstable. 



For Celn3, Lawrence et al. [jl6[ interpreted the deviation of x(T) from the Curie- Weiss 
behavior in the low temperature region as arising from spin-fluctuations with the character- 
istic temperature Tgf ~ 50 K. We believe that this effect is clearly from the CEF splitting 
(A78=125 K). The calculated Tx (few K) is much lower than CEF splitting and is inacces- 
sible experimentally due to the onset of anti-ferromagnetism at 10 K. This CEF effect has 
been checked by comparing the experimental results with the Curie- Weiss susceptibility of 
Fy- and Fg— magnetic moments as well as the van Vleck susceptibilities. In particular, the 
experimental low temperature x{T) curve shows saturates only below 10 K Tgf) and 
above this is well-fitted to the Curie- Weiss form with a characteristic temperature 18 K, 
which can be interpreted as a combination of Neel temperature plus the CEF effect from 
the excited quartet field Fg. 

In FIG. 3(b), susceptibility curves of CeSna are plotted for several values of e/ with a 
comparison to experiment [|T8|. Since xiT 0)(~ ^/Tk) is exponentially sensitive to e/ 
(see Eq. (4)), we performed calculations for several e/'s (—2.0, —2.1, —2.2 eV). Clearly our 
calculations bracket but do not fit the experimental data, since a calculation for A78=0 
fits the data well we suspect the source of the disagreement may be an overestimate 
of crystal field splitting in our calculation, placing the effective degeneracy between two 
(Tk < Ays) and six (Ayg = 0). 

In conclusion, we reproduced both the band and many body features of the Ce- 
compounds by inputting LDA(LMTO-ASA) calculated hybridizations into f/ — oo NCA 
calculations for appropriately defined impurity Anderson models. We reproduce well exper- 
imental trends in Tk-, Ayg, and x{T) for the CeMa series. This work provides a starting 
point for quantitative calculations of realistic heavy fermion systems at finite temperatures. 
Improvement may come through a proper inclusion of dynamics [U ^ oo), a reliable 
theory for electrostatic CEF contributions, and through self-consistent closure of the many 
body calculations. 

We are grateful to J. W. Allen, M. Steiner and J. W. Wilkins for helpful discussions. 
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FIGURES 

FIG. 1. Schematic diagrams of hopping of /-electrons to band states in Cea;Lai_2;M3 systems, 
(a) In the X = 1 hmit, the center /-orbital couples the rest of the lattice (enclosed in dashed line), 
i.e., effective medium. The /-orbital is connected to ligand orbitals (cj) with Vim (i lattice site, m 
the CEF index) . The lattice of / orbitals (origin excluded) are treated as producing a static mean 
field coupling to ligands in the LDA calculation, (b) In the impurity limit (x ^ 0), the /-orbital 
(at center) couples to the bare ligand bands. 



FIG. 2. Hybridizations of Cea;Lai_a;M3 calculated in two different limits of x. (a) x = 1 limit. 
Although Tt{J = 5/2) (thick line) and Ts{J = 5/2) (thin line)-hybridizations are comparable, large 
hybridization (/^ <r^ fluctuation) for Fy near 1.0-3.0 eV pushes the doublet below the quartet. 
The spikes appearing above the Fermi energy (dashed line) are due to the coupling of /-orbital at 
origin to the bonding (anti-bonding) states of ligand and neighboring /-orbitals. (b) a; — > limit. 
Curves are almost identical to (a) except near the Fermi energy. Note that the peak at —0.1 eV 
(inset) is pushed down to —0.2 eV (inset (a)) due to /-ligand coupling. 



FIG. 3. Calculated magnetic susceptibility x{T)-vs.-T for CeM^ compounds, (a) M dependence 
of x(T); see text for discussion, (b) Magnetic susceptibilities of CeSna for different /-level positions. 
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TABLES 

TABLE I. Table of the electronic configuration of ligand atom M (E.G.) ground CEF multiplet 

Tgrd, CEF splittings A-jg and Kondo temperatures Tk for Cca;Lai_a;M3 systems. Units are in 
Kelvin. Positive A78 indicates a stable ground doublet on the Ce sites. (*: For mixed valence 

systems, CEF's are not experimentally dedueibl(^ siuee Ay^s -C -/a)- 
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(a) concentrated limit (CeM3) 




(b) impurity limit (Ce La._M„,x^O) 




-12 -8 -4 4 8 

E eV 



10 



1-1 



o 
E 



10 



1-2 



E 



— I — I — I 1 1 1 1 1 1 — I — I — I 1 1 1 1 1 1 — I — I — I 1 1 1 1 1 1 — I — I — I 1 1 1 1. 

A : 



^ 



^ 



^ A 



Curie— Weiss 



■(a) 



o CeSn3 
A CePb3 
Celn3 



OOOGOOOO 



O O Q 



■ ■■■■■■ I I ■■ l_ 



■ 



0.1 



10 



100 



1000 



10000 



3.0 







O 




E 


2.4 


=5 




E 






1.8 








1.2 






O 





0.6 



CeSn3 





• Experiment 
o Theoryref=-2.0eV) 
A Theory(£;f=-2.1eV) 
Theory(ef=-2.2eV) 



• 

A 



(b) 



200 



400 

T (K) 



600 



800 



